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ABSTRACT 


Finite element analysis is applied to study the large deflection and 
stress of a thin-walled pressurized torus loaded hy normal contact with a 
plane. The torus is found to have an elliptical footprint area, and con- 
siderable bulge occurs in the sidewall in the vicinity of the load plane. 
In large load ranges, the finite element calculations show compressive 
circumferential stress and negative curvature in the footprint region. 

An experimental study of the standing torus, using liquid metal strain 
gages , is outlined. Experimentally determined stresses are compared to 
those resulting from finite element analysis at various meridional and 
circumferential coordinates of the torus , including the footprint area. 
Circumferential strains compare favorably while meridional strains are 
higher in the finite element analysis, probably due to slippage of the 
boundary at the rim. 

This study utilized the STAGS finite element computer program. The 
purpose of the study was to evaluate the various program options for 
structural loading, for material modeling, for stress analysis, and for 
grid refinement. The experimental model, a thin-walled rubber tube mounted 
on a steel cylindrical rim, provided measured results to compare with 
various analytical trade-offs. It was found that there was almost no 
difference in predicted deflections or stress distributions between linear 
and nonlinear material description. However, the difference between lin- 
ear analysis and that of geometric nonlinearities utilizing incremental 
loading was marked . 


TORUS DIMENSIONS AND GLOBAL COORDINATE SET 


STAGS (ref. l) uses a first order shell theory to reduce a 3- 
dimensional structural problem to a dependency on 2 spatial coor- 
dinates , here the meridional coordinate a and circumferential 
coordinate 9 . The deflection of the plane into the torus is 

5 . An example of an 11 x 15 torus grid is shown and dimen- 

sions of the cylindrical rim and inflated rubber tube are given. 
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EXPERIMENTAL SET-UP SHOWING ELLIPTICAL FOOTPRINT 

The tube was inflated to 0.03^ atm. (0.5 psi) and test loads of 
3k. 3 N (7.7 Us), 69.9 N (15.7 Its), and IlL.L N (25.7 Its) were applied. 
Deflection, footprint area, circumferential and meridional strain, and 
final pressure were ot served. 
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LIQUID METAL STRAIN GAGE AND CALIBRATION DEVICE 


In order to measure the large strain encountered , strain gages with 
a rubber housing and capillary column of mercury were used. The gages 
are inherently nonlinear by nature, and individual calibration is re- 
quired (ref. 2). Derivation of gage factor versus strain is given in the 




GAGE PLACEMENT AT 0° MERIDIONAL ANGLE - 
COMPLETE SYSTEM WITH DATA AQUISITION 



Calibrated gages were mounted on the torus with silicone rubber 
cement. The torus was inflated to test pressure and rotated to the de- 
sired circumferential angle, load was applied, and readings were re- 
corded. The microprocessor shown below tabulated and recorded, on 
magnetic tape, data for strain versus circumferential angle at given 
values of load and meridional angle. 
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DATA ACQUISITION FLOW CHART 


The logic flow chart for the data acquisition scheme is shown "below. 
The buffer for the strain signals and the program controller with the 
stored logic of the acquisition system are encompassed in the micro- 
processor. The microprocessor is programmed in the BASIC language, and 
the digitizing of analog gage voltage is accomplished by a digital volt- 
meter. The three channels referred to are the meridional strain, the 
circumferential strain, and a dummy gage voltage. 
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LOAD-DEFLECTION CURVES FOR ANALYTICAL FINITE ELEMENT 

AND EXPERIMENTAL DATA FOR TORUS AT P. = 0.0 3^ ATM 

l 

This figure shows points on the load-deflection curve obtained from 
a STAGS analysis using a 13 x 15 quarter torus at loads of 3^.3 N 
(7.7 lbs), 69.9 N (15.7 lbs), and llU.L N (25-7 lbs). The deflections 
predicted from STAGS are lower since minimum potential theory is em- 
ployed and structural stiffness is overestimated with a coarse grid. 

As grid size is refined, predicted deflections approach the experi- 
mental values, but computer run time also increases (ref. 3). At a 
load of 3^.3 N, three different grid sizes were modeled. 
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' 13 x ; 15 QUARTER TORUS GRIDS WITH 
■ 3^.3 N (.7.7 Tbs ) LOAD - NONLINEAR ANALYSIS 

This figure gives the computer-generated plot of the quarter torus with 
a 13 x 15 grid. Exploitation of the torus symmetry allows analysis of 
one-fourth of the structure, giving a grid density four times greater than 
that of a 13 x 15 grid applied to the full torus. The unloaded grid geom- 
etry is overplotted . on the the deformed geometry. Note the plane surface 
at the load point (footprint) and the side wall bulge. The nonlinear 
geometric analysis option of STAGS gave convergence to a final deflection 
of 8k% of the experimental value while the linear analysis yielded only 
21# of the experimental value. This result is significant when compared 
to toroidal analysis using other finite element procedures (refs, k and 5). 




CALCULATED STRESS DISTRIBUTION AT 90° MERIDIONAL 
ANGLE FOR 13 x 15 QUARTER TORUS GRID - L = 3^.3 N (7*7 lbs.) 

Calculated values of meridional stress and circumferential stress 
versus circumferential angle 0 are shown below. The variation of 
stress in the footprint and bulge regions is significant. The values 
remain constant away from the contact area and are mainly due to in- 
ternal pressure. These values, with this grid refinement, agree 
favorably with other investigators (ref. 6). 





CALCULATED MERIDIONAL STRESS DISTRIBUTION, ° a , 13 x 15 

QUARTER TORUS GRID - NONLINEAR ANALYSIS, LOAD = 34.3 N (T - T IDs) 

The three dimensional plot below shows meridional stress dis- 
tribution over the surface of the torus at a 34.3 N load. Even 
tensile stress appears in both the meridional and circumferential 
directions in the upper half of the torus , as a result of internal 
gas pressure. However, the compressive effect of contact with the 
flat plate reduces these tensile stresses in the footprint region. 
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CALCULATED CIRCUMFERENTIAL STRESS DISTRIBUTION, °0 , 13x15, 
QUARTER TORUS GRID - NONLINEAR ANALYSIS, LOAD = 3U.3 N (j.j lbs) 

A three dimensional plot of circumferential stress distribution 
over the upper half of the torus is shown below. As with the previous 
figure, uniform tensile stress due to gas pressure is significantly 
reduced in the footprint region, to almost 0 N/mm 2 in the center. Results 
of an intermediate load step of 6 9-9 N (15.7 Lbs) are contained in 
the literature (ref. 7). 
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13 x 15 QUARTER TORUS GRID WITH 
11U.U N (25.7 lbs) LOAD - NONLINEAR ANALYSIS 

The figure below shows the STAGS generated plot of a quarter 
torus in the unloaded and deformed configurations under a test 
load of 11^.^ N. The deflection has increased significantly and 
the degree of bulge in the sidewall is more pronounced. The center 
of the footprint region shows negative curvature. 


CALCULATED CIRCUMFERENTIAL STRESS 
DISTRIBUTION, a e , AT LOAD = 1 1 U . 4 N (25-7 Its.) 

A three dimensional plot of circumferential distribution over 
the quarter torus is given belov for the llU.U N load. Note the 
scale change from that of the previous two three dimensional plots . 
Increased compressive stress in the footprint region is seen to cause 
a stress reversal in the center of the footprint region. The stress 
becomes compressive near the center, then suddenly reverses and 
becomes tensile, suggesting that a local limit point on the structure 
load deflection path has been reached and snap-through has occurred. 


3.5 x 1CT 2 N/mm 2 (50 PSD 



FOOTPRINT REGION 



FOOTPRINT OF AN INNER TUBE SHOWING NEGATIVE CURVATURE 


n* ,Gf HAL * 


The figure below is a photograph of an inner tube under rela- 
tively heavy load. An initial eccentricity has been introduced. 
The negative curvature near the center of the footprint is consid- 
erably reduced in the presence of friction at the contacting sur- 
faces (not modeled here). 
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Strain (percent) 


EXPERIMENTAL AND CALCULATED STRAINS AT 
-10° MERIDIONAL ANGLE - LOAD = 3^.3 N ( T - 7 lbs.) 

The figure below- compares measured and calculated meridional 
and circumferential strains for varying 0 at « = -10°. This -10° 
meridional angle places the strain gages in the region near the rim. 
Good agreement is achieved in observed and calculated values of 
circumferential strain. The disagreement between observed and cal- 
culated meridional strain is attributed to slipping of the membrane 
under the rim on the test fixture. The computer model had a fixed 
boundary at the rim. 


5 Meridional Strain (Predicted) 

(£) Meridional Strain (Measured) 

6 Circumferential Strain (Predicted) 



0 10 20 30 40 50 60 70 80 90 


Circumferential Angle (Degrees) 


137 



EXPERIMENTAL AND CALCULATED STRAINS AT 
90° MERIDIONAL ANGLE - LOAD = 69-9 N ( 1 5 - T lbs.) 

This figure compares measured and calculated strains for vary- 
ing circumferential angle 0 at a = 90°. At this farthest point 
from the rim, meridional strains are in good agreement. Note here 
that for all measurements below a = 30°, the strain gage is in the 
footprint region. Measured strain did not become compressive here, 
probably due to friction between the contacting surfaces. Reports of 
strain comparison at other loads are contained in the literature 
(ref. 8) . 
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APPENDIX 


In describing the operation of a liquid metal strain gage, one 
assumes that the liquid metal column obeys the fundamental resist- 
ance relation (there are no voids in the column). Further, the 
resistivity of the conducting medium is assumed constant, and the 
volume of the capillary cylindrical cavity is assumed constant. 
Rubber is essentially incompressible; the assumption of constant 
volume under elongation is a valid one. The fundamental resistance 
is 
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when p = resistivity, ohm-in. 

L 2 length, in. 

A = cross-sectional area, in. 
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Using a Taylor series expansion about the initial length, 
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where e = strain, in. /in. 

Equation 2 establishes a quadratic relationship between resist- 
ance change and strain. It is seen from ( 2 ) that for small strains, 
the gage factor is two. 
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